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I.  INTRODUCTION 


The  traditional  structure  used  in  containing  very 
high  pressures  is  the  thick-walled  cylinder.  The  solution 
for  the  stress  in  these  vessels  was  developed  by  Lame  and 
is  documented  in  most  primary  texts  in  Strength  of  Materials* 

I 

The  Lame  solution  reveals  that  the  highest  tensile  stressed 
area  in  thick-walled  cylinders  is  the  inside  (bore)  radius. 

To  reduce  the  bore  stress,  several  techniques  have  been 
developed  to  induce  compressive  residual  stresses  at  the 
bore  and  thus  increase  the  pressure  that  can  be  withstood 
elastically.  They  are  multi-layer  construction,  autofrettage 
and  wrapping2.  The  primary  concern  of  this  report  is  the 
autofrettage  method. 

Autofrettage  is  a process  where  a thick  walled  cylinder 
is  subjected  to  a pressure  sufficient  to  cause  yielding  in 
the  cylinder.  The  pressure  is  then  released  and  the  cylinder 
recovers.  Since  the  outside  portion  of  the  cylinder  was 
subjected  to  small  amounts  of  strain,  it  will  try  to  return 
to  nearly  its  original  dimension,  while  the  bore  is  subjected 
to  relatively  large  plastic  strain,  and  will  try  to  remain 
deformed.  This  non-uniform  recovery  causes  a compressive 
residual  stress  at  the  bore  and  a tensile  residual  stress  at 
the  outside  radius.  Figure  1 is  a plot  of  some  typical 
autofrettage  residual  stress  distributions. 
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TYPICAL  AUTOFRETTAGE  RESIDUAL 
STRESS  DISTRIBUTIONS 


^ 100%  OVERSTRAIN 


BORE  RADIUS  OUTSIDE  RADIUS 


FIGURE  I. 


The  maximum  compressive  residual  bore  stress  is  obtained 
when  the  autofrettage  pressure  is  sufficient  for  the  outside 
radius  to  yield.  This  is  called  the  100%  overstrain  con- 
dition*. This  condition  allows  for  the  maximum  pressure 
to  be  contained  elastically  and  has  traditionally  been  con- 
sidered the  optimum  design  condition.  If  an  autofrettaged 
cylinder  is  subjected  to  fatigue,  the  optimum  autofrettage 
condition  may  not  be  100%  overstrain.  The  combination  of 
tensile  residual  stress  and  tensile  operating  stress  at  the 
outside  radius  may  cause  crack  initiation  at  this  surface. 
This  report  is  the  result  of  an  investigation  of  the  effects 
of  the  autofrettage  residual  stress  distribution  on  crack 
propagation  initiating  from  flaws  at  the  outside  radius. 

I I . APPROACH 

In  order  to  model  the  effect  of  autofrettage  on  external 
flaw  initiated  fatigue  crack  propagation,  a model  must  first 
be  developed  for  the  crack  propagation  in  a cylinder  that 
has  not  been  autofrettaged.  This  can  be  accomplished  with 
the  use  of  linear  elastic  fracture  mechanics,  and  the  Paris 
power  law  for  crack  propagation^: 

^ = CjAkJ  (la) 

dN 

‘Percent  overstrain  is  defined  as  the  percentage  of  the 
wall  thickness  that  is  subjected  to  plastic  strain  during 
auto  f re  1 1 age . 


Since  the  primary  purpose  for  modeling  crack  growth 
is  to  determine  remaining  fatigue  life,  equation  la  can  be 


integrated  to: 

cf 

(Nf-Ni)  = J _dc (lb) 

c.  C ! A Km 

It  is  assumed  that  autofrettage  will  have  an  effect  on 
crack  propagation.  This  effect  may  be  modeled  by  the  following 
modification  to  the  power  law: 


dc 

dN 


H'CjAK11’, 


(2a) 


where  ¥ is  an  arbitrary  function  of  autofrettage  parameters 
determined  from  experiment.  Once  4*  is  known,  an  expression 
for  the  remaining  fatigue  life  of  an  autofrettaged  cylinder 
can  be  determined  by  integrating  equation  2a: 


(Nf-Ni) 


/ 


dc 


fC  j A Km 


(2b) 


Based  on  the  above  rationale,  it  was  planned  to  model 
the  crack  propagation  by  determining  stress  intensity  factors 
at  various  crack  depths  in  non-autofrettaged , externally  flawed 
thick-walled  cylinders  subjected  to  internal  pressure.  It 


4 


was  also  planned  to  determine  the  function  V by  performing 
the  simple  experiment  discussed  below. 

III.  THEORY 

Several  investigators  have  attempted  to  determine  the 
stress  intensity  factors  for  an  externally  flawed  cylinder 
subjected  to  internal  pressure.  Kobayashi4,  Chopra^,  and 
Kobayashi  et  al6,  have  found  the  stress  intensity  factors  in 
thin-walled  cylinders  using  various  techniques.  Emery  and 
Segedin7  have  determined  the  stress  intensity  factors  for 
cylinders  of  wall  thickness  to  inside  radius  ratio  (B/a), 
varying  from  .052  to  1.85  using  a finite  difference  method. 
While  this  method  yields  good  results  for  some  loading  and 
specimen  geometries,  it  gives  poor  results  in  others.  It  was 
therefore  decided  to  develop  an  expression  for  the  stress 
intensity  factor  using  another  method. 

One  assumption  that  can  be  made  is  that  a single  external 
crack  in  a thick-walled  cylinder  will  produce  the  same  stress 
intensity  factor  as  that  encountered  in  a single  edge 
notched  panel.  The  stress  intensity  factor  can  then  be 

O 

approximated  by  : 

K = 1 .12a/  c tt  , (3) 


where  a is  the  stress  normal  to  the  piano  crack,  in  the  vicinity 
of  the  crack. 
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The  problem  then  reduces  to  evaluating  the  tangential 
stress  in  the  uncracked  ligament  of  wall  thickness.  This 
stress  can  be  approximated  based  on  the  following  assumptions 

1.  The  cylinder  is  infinitely  long  and  open-ended. 

2.  The  crack  is  oriented  in  the  radial  direction  and 
is  of  constant  depth  for  the  entire  length  of  cylinder. 

3.  The  tangential  stress  distribution  (at°°)  , 180°  re- 
moved from  the  crack  is  unaffected  by  the  presence  of  the 
crack  and  is  given  by  the  Lame  solution*. 


a 


too 


[1  - -2] 
T-  2 


(4) 


4.  The  tangential  stress  distribution  in  the  uucracked 
ligament  of  wall  thickness  (°tmod^  *s  given  by  the  Lame 
solution  for  an  uncracked  cylinder  whose  outside  radius  has 
been  reduced  by  the  crack  depth: 


atmod 


P ; a ‘ 


(b-c)^  - a ^ 


[ * 


+ 


(b-c)2 


(5) 


5.  The  stresses  are  symmetric  about  the  plane 
The  first  two  assumptions  reduce  the  analysis  to 


of  the  crack 
a two- 


dimensional  problem.  The  third  is  an  evocation  of  the  St. 


Venant  Principle.  The  fourth  is  a "uniform  thinning" 
approximation  that  has  been  used  successfully  in  analyzing 
internally  cracked  cylinders9,  and  the  last  is  based  on 
the  geometry  of  a cracked  cylinder.  These  assumptions  are 
represented  schematically  in  Figure  2a. 

The  stability  of  a cylinder  subjected  to  the  above  assumed 
loading  can  then  be  determined  by  analyzing  its  free  body 
diagram.  Figure  2b.  It  is  readily  observed  that  the  cylinder 
is  in  translational  equilibrium.  Since  the  resultant  forces 
of  the  stress  distributions  act  eccentric  to  the  resultant 
force  of  the  pressure  distribution,  a moment  (Mg)  must  act 
over  the  uncracked  thickness  for  rotational  equilibrium  to  be 
maintained.  The  value  of  Mg  per  unit  length  of  cylinder  can 
be  shown  to  be: 


(b-c)2  b-c 

— In  [ ] 

(b-c) 2-a2  a 


(6) 


The  total  tangential  stress  distribution  in  the  uncracked 


wall  thickness  therefore,  has  two  components:  1)  the  modified 
Lame  stress  and  2)  a bending  stress  due  to  Mg  (ot  bend) • The 
bending  stress  can  be  found  based  on*: 


MBy 

°t  bend  = ’ (7a) 


7 


Figure  2a. 


SI^’SSSS  °o 


Figure  2b.  Free  body  diagram  of  an  externally  flawed,  thick-walled 
cylinder,  subjected  to  internal  pressure 


where  y is  the  distance  from  the  neutral  axis  and  I is  the 
moment  of  inertia.  Since  the  area  over  which  Mfi  acts  is 
rectangular  as  shown  in  Figure  2c,  the  bending  stress  in 
this  case  is: 


1 2P ■ a2 [r-  (h-c  + a) j 
°t  bend  = 2 ' k 


(b- c- a) ^ 


In  [-] 


<b-c)2  1„[^£ 


.b2-a2  a (b-c)2-a2  a 

(7b) 

The  total  tangential  stress  distribution  is  therefore 
given  by  the  sum  of  equations  7b  and  5: 


Of  * P 


i»2/ - n + (b-c)  2]f 

\(b-c) 2-a2  r2 


12  (r‘ 


(b- c+a) 


(b-c-a)3  Lb2-a2 


) T b2  fbl 

- — In  - 

Lb2-a2  La. 


(b-c) 


(b-c) 2-a2 


In 


(8) 


To  check  the  accuracy  of  equation  8,  it  was  decided  to 
utilize  the  finite  element  program  NASTRAN.  A grid  was 
generated,  using  linear  elements,  to  model  a cylinder  of  4 
units  outside  radius  and  2 units  inside  radius,  for  a B/a 
ratio  of  1.  The  grid  is  illustrated  in  Figure  3.  Cracks  of 
various  depths  were  modeled  by  running  the  program  several 
times.  On  each  successive  run,  the  constraint  on  the  row 
of  grid  points  marked  A-A  was  reduced,  which  allowed  the 
unrestrained  points  to  act  as  a free  surface.  The  NASTRAN 


1 


m 


predicted  stresses  are  plotted  in  Figure  4 and  show  excellent 
agreement  with  equation  8. 

The  stress  intensity  factor  is  determined  by  sub- 
stituting the  quantity  (b-c)  for  r in  equation  8,  and  using 
this  value  of  a t for  0 in  equation  3: 


K = 1 . 12P  . a^/iTc 
I i 


. 2 2 
(b-c)  -a 


6 (b- 
(b  - c 


c + ajT 

-a)n, 


•b2-a2 


1 n [— ] 
a 


(b-c) 2-a2 


(9) 


The  range  of  stress  intensity  factor  during  a fatigue 
cycle  is: 

A K i = 1 . 1 2 A P i a 2 / ir  c 


[— - 
L (b-c)  -a  (h-c-a)5^  -a^ 


(b‘c)  - 


2 2 

(b-c)  -a 


(10) 


This  expression  is  substituted  in  equation  lb,  to 

yield  the  following  integral  for  remaining  fatigue  life 
_cf 


(N 


f-N.)  = c^^l.12  AP.a2  /Vc 
ci 

[2  6 (b- c+a)  T b2 

(b-c) ^-a2  (b-c-a)^b2-a2 


b<-  , rbT 

ln[g] 


^ 2 2 
(b-c)  -a 


(ID 


Since  equation  11  can  only  be  integrated  by  approximation, 
it  was  decided  to  develop  an  analytical  expression  for  re- 
maining fatigue  life,  by  determining  the  stress  intensity 
factors  by  another  method.  A recent  study  by  Hussain  et  al^ 
has  found  that  stress  intensity  factors  can  be  predicted  very 
accurately  and  efficiently  by  using  quadratic  finite  elements 
in  the  NASTRAN  program. 

The  element  is  basically  a trapezoid  that  is  defined  as 
the  connection  of  eight  grid  points,  (four  corner  points 
and  four  mid-side  points)  as  shown  in  Figure  5a.  If  one  side 
of  the  trapezoid  is  defined  as  connecting  two  corner  and  one 
mid-side  points,  all  having  the  same  coordinate  position, 
these  points  define  one  corner  of  a triangle  as  shown  in 
Figure  5b.  If  the  adjacent  mid-side  points  are  then  placed 
one  fourth  of  the  distance  between  the  remaining  corners  of 
the  triangle,  the  numerical  expressions  for  the  stresses  and 
strains  at  these  mid-side  points  are  singular  of  order  one 
half  (l//r),  which  is  the  same  order  of  singularity  that 
occurs  in  the  Westergaard  near  field  solution11.  These 
elements  can  therefore  be  used  to  very  accurately  model  a 
crack.  The  stresses  or  strains  predicted  by  NASTRAN,  using 
these  elements,  are  substituted  in  the  Westergaard  solution 
to  solve  for  the  stress  intensity  factor.  This  method  has 
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FIGURE  6 QUADRATIC  FINITE  ELEMENT  GRID 


I.« 
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predicted  stress  intensity  factors  for  several  geometries 
within  2-3%  . Figure  6 is  a grid  with  six  singular 
elements . 

Stress  intensity  factors  were  found  for  nine  crack 
depths.  The  results  are  compared  with  equation  9 in 
Figure  7.  The  plot  shows  that  equation  9 predicts  stress 
intensity  factors  much  greater  than  the  NASTRAN  prediction 
for  deep  cracks.  A deviation  from  equation  9 can  be  ex- 
pected for  deep  cracks,  since  it  was  derived  from  equation 
2,  which  was  developed  assuming  that  the  crack  depth  is 
negligible  in  comparison  with  the  size  of  the  medium  in  which  it  is 
embedded®.  For  very  shallow  cracks,  the  assumption  is  valid, 
but  for  deep  cracks  it  does  not  hold;  therefore,  the  NASTRAN 
predictions  are  believed  to  be  far  more  accurate  than  equation  9. 

To  fully  utilize  the  NASTRAN  data,  it  was  decided  to 
fit  a curve  to  them.  An  equation  of  the  following  type  was 
assumed: 

KI  = 0 ^ Y[|]  . (12a) 

where  a is  a stress  far  removed  from  the  crack. 

For  the  case  of  internal  pressure  loading  o was  chosen 
as  the  tangential  stress  at  the  outside  radius  180°  removed 
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from  the  crack  (a  = °too  @(r=b)).  The  expression  for  the 
stress  intensity  factor  becomes: 


2Pia* 

T^I2 


/B  Y r- 


[-] 

B 


(12b) 


The  arbitrary  function  Y [£]  is  found  by  linear  re- 
gression. The  best  fit  is  accomplished  by  the  fourth  de- 
gree polynomial: 

Y[l]  = .5687  + 18. 3132 [“] 2 - 32.7075[^]3  + 25.5121[|]4 

(13a) 

2 

This  function  has  a correlation  coefficient  (R  ) of  .9999. 

When  equations  13a  and  12b  are  substituted  in  equation  la 
to  predict  remaining  fatigue  life,  the  resulting  expression 
is  very  cumbersome  to  integrate. 

A good  fit  can  also  be  obtained  by  the  quadratic: 

Y [— J = .6916  + 8.1217[|]2  (13b) 

B D 

which  has  a correlation  coefficient  of  .9911.* 

Substituting  equation  13b  into  equation  12b  the  following 
is  obtained: 

2 P • a 2 c 

Kt  = 1 /B  (.6916  + 8.1217 [— ] 2 ) (14) 

b 2-  a 2 B 
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The  range  of  stress  intensity  factor  during  a fatigue 


cycle  is  given  by: 


2A  P . a 2 _ 

AK  = /B  (.6916  ♦ 8.1217  2 ) 

1 b2-a2  B 

(15) 


When  equation  15  is  substituted  in  equation  lb,  the  ex- 
pression can  be  integrated  in  closed  form  for  integer  values 
of  Paris  power  law  exponents  m.  Since  most  steels  have  m 
values  between  2 and  4,  the  integration  has  been  performed 
for m = 2,3  and  4.  The  results  of  the  integration  are  re- 
ported in  Table  1.  Equations  13a  and  13b  are  compared  with 
the  NASTRAN  data  in  Figure  8a. 

The  NASTRAN  predicted  stress  intensity  factors  are 

compared  with  those  predicted  by  Emery  and  Segedin^  in 

Figure  8b.  The  plot  shows  that  the  NASTRAN  prediction  differs 

substantially  from  that  predicted  by  Emery  and  Segedin  for 

thick-walled  cylinders  (B/a  >.5),  but  is  still  in  trend 

agreement  with  their  predictions  for  thinner  walled  cylinders 

(B/a  = .111,  .25).  There  is  no  reason  to  believe  the  NASTRAN 

data  to  be  in  error,  and  it  is  therefore  assumed  to  be  the 

best  solution  available  for  cylinders  of  B/a  ratio  of  1. 

IV.  EXPERIMENTATION 
— 

To  determine  the  effects  of  the  autofrettage  residual 

. 


I 


let 


m = 


ra 


TABLE  1.  Expressions  for  remaining  fatigue  life  in  an  externally 
flawed,  thick -walled  cylinder,  subjected  to  internal 
pressure  for  various  values  of  Paris  power  law  exponent 


(Nf-Ni)  = 


cf 


2APia‘ 


(b2.a2)B3/2 


r 


dc 


6916B2  + 8.1217c2) 


m 


ci 


$ = 


U 


2 Pia^ 
(b2-a2)B 


a = .6916B2 

6 = 8.1217 
2 


(Nf-N 


i ) = $ j” 


2a(a+6c2)  2/a3 


1 f 

1 t a n " * c /a  6 S 

a3  a **  ci 


(Nf-Ni) 


rc  (5a  + 3 3c2)  + 3 tan'  zJ 

L8oi2  (a+Bc272  8 /a ^8 


/a  3 
a 


ci 


= 4 


cf 


(Nf-Nj)  = * 


rC(33a2'l~40a3c2  + 15  82c4)  + 15  tan'1  c/ag  j 

L 48a3(a+Bc2)  3 48^T8  a 

ci 


i 
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NORMALIZED  CRACK  DEPTH  (c/B) 


FIGURE  8o  NASTRAN  STRESS  INTENSITY  FACTORS 


■ •>  . •> 

l-.r.  - — WJ<r. < 


normalized  crack  OEPTH  (c/B) 


FIGURE  8b  COMPARISON  OF  STRESS  INTENSITY  FACTORS 
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stresses  on  crack  propagation  and  thus  the  function  y,  a 
simple  experiment  was  proposed  using  thin  disks  of  miniature 
autofrettaged  cylinders  . The  disks  were  loaded  diametri- 
cally and  had  machined  notches  at  the  outside  radius  as  shown 
in  Figure  9.  Three  amounts  of  autofrettage  were  tested;  0%, 
50%  and  100%  overstrain.  The  material  that  the  disks  were 
made  from  had  the  following  chemistry  and  properties  (AISI 
4340)  : 


% by  weight 


C Mn  Si  Ni  Cr  Mo  . Fe 

.39  .80  .29  1.97  .8  .25  Bal. 


.1%  Yield  strength  166,000  psi  (1144  MPa) 

Ultimate  strength  175,000  psi  (1206  MPa) 

The  specimen  dimensions  were: 

a = . 604  in . (1.53  cm) 

b = 1.125  in.  (2.86  cm) 
t = .25  in.  (.635  cm) 

The  applied  alternating  load  (AF)  was  2700  lb.  (12000N). 


In  the  autofrettaged  specimens  (50%,  100%  overstrain) 
fatigue  cracks  readily  initiated  at  the  root  of  notch.  In 
order  to  initiate  cracks  in  the  non- autofrettaged  specimens 
(0%  overstrain),  much  deeper  notches  were  required. 

V.  RESULTS  AND  DISCUSSION 


To  compare  the  results  of  the  experiment  on  a c vs,  N 
type  plot,  N had  to  be  arbitrarily  defined  as  the  number 
of  applied  cycles  minus  the  number  of  cycles  required  to 


grow  a crack  1/2  of  the  distance  through  the  wall  thickness 
(N-N50) . The  raw  c vs.  (N-Njq)  data  is  shown  in  Figure  10. 

An  interesting  phenomenon  occurs  for  very  deep  cracks. 

The  crack  growth  rate  slows  in  the  radial  direction.  and  the 
crack  gradually  changes  direction  and  grows  parallel  to 
the  applied  alternating  load.  This  phenomenon  occurs  when 
a crack  is  first  subjected  to  Mode  I and  then  Mode  II  defor- 
mations, and  has  been  investigated  by  Hussain  et  al*3. 

The  loading  of  the  disk  as  described,  produces  a compressive 
normal  force,  and  a moment  on  the  plane  ahead  of  the  crack. 

The  moment  causes  a tensile  bending  stress  sufficiently  large 
such  that  the  resultant  of  the  bending  and  compressive  normal 
stresses  is  still  tensile  in  the  outer  members  of  the  disk, 
in  the  direction  perpendicular  to  the  crack.  This  tensile 
stress  causes  the  Mode  I opening.  As  the  crack  grows,  the 
line  of  action  of  the  normal  force  becomes  closer  to  the 
line  of  action  of  the  applied  force,  the  moment  required 
for  equilibrium  therefore  becomes  smaller.  For  very  deep 
cracks,  the  normal  stress  becomes  dominant  causing  Mode  I 
closure.  The  compressive  stress  then  causes  a shearing  of 
the  crack  tip  or  Mode  II  sliding,  which  causes  the  crack  to 
change  direction. 

To  determine  T , the  stress  intensity  factors  had  to  be 
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determined  for  non- autof rettaged  disks  under  this  loading 
geometry.  This  was  accomplished  by  again  using  NASTRAN 
with  singular  elements.  The  same  grids  were  used  as  with 
the  internal  pressure  loading.  To  present  the  results,  a 
relation  of  the  form  of  equation  12a  was  again  assumed. 

In  this  case,  the  stress  a was  taken  to  be  the  absolute 
value  of  twice  the  normal  stress  (a  = F/Bt),  therefore: 

K = _L_  /B  Y [— ] (16) 

1 Bt  B 

The  best  fit  to  the  NASTRAN  data  revealed: 

Y[|]  = . 3625+2. 1834 [|]-2. 8386 [£] 2- 1 . 0438[|] 4 

Equation  17  is  compared  with  the  NASTRAN  results  in 
Figure  11.  This  plot  predicts  the  slow  crack  growth  for  deep 
cracks,  because  of  decreasing  stress  intensity  factor. 

Equations  16  and  17  were  used  to  develop  a log  dc/dN  vs. 
log  AKj  plot  for  the  non-autofrettaged  disks.  The  expected 
straight  line  behavior  was  not  observed  due  to  the  small 
range  of  A Kj  developed  in  this  test  (A  Kj < 10  KSI*in.). 

The  straight  line  correlation  occurs  at  higher  values  ofA  Kj 
(20  KS  l **1  n . < A Kj  < 80  KSI  /in  . ) , and  is  often  extrapolated 
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through  the  lower  range.  The  data  does  fall  within  the 
scatter  band  that  has  been  developed  for  434014  and  is 
shown  in  Figure  12.  This  experiment  therefore,  could  not 
be  used  to  determine  the  function  V . 

The  raw  c vs.  (N - N 5 Q ) data  (Figures  13,14,15)  do  yield 
interesting  results  however.  When  the  curves  for  the  three 
amounts  of  overstrain  are  superimposed  (Figure  16),  it  can 
be  readily  seen  that  autofrettage  definitely  increases  crack 
growth  rate  in  the  disk  specimens.  The  magnitude  of  this 
effect  can  be  better  seen  in  Figure  17  which  plots  crack 
growth  rate  (dc/dN)  vs.  crack  depth.  The  plot  shows  crack 
growth  rate  is  substantially  increased  with  50%  overstrain 
and  further  increased,  but  to  a lesser  extent,  with  100% 
overstrain.  This  drastic  change  was  not  expected. 

The  increase  in  crack  growth  rate  may  be  due  to  the 

tensile  residual  stresses  in  the  outer  members  of  the 

cylinder  increasing  the  mean  stress  rra  +0  -i  / 2)  and 

min  maxJ 

the  stress  ratio  (0  /o  )•  It  is  well  known  that  crack 

min  max 

growth  rate  increases  with  increasing  mean  stress  and  stress 
ratio,  but  generally  not  to  the  extent  that  was  observed  in 
these  experiments. 


In  examining  the  specimens,  it  was  noticed  that  in  the 
non  - auto f re tt aged  specimens,  crack  growth  was  not  accompanied 


by  visual  plastic  deformation,  while  in  both  the  50%  and  100% 


FIGURE  14  c VS.  N (50%  OVERSTRAIN) 
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overstrain  specimens,  plastic  zones  were  readily  visible. 

The  plastic  zones  were  noticed  when  the  crack  was  about 
38%  through  the  wall  thickness  in  the  50%  overstrain 
specimens  and  when  the'crack  was  about  31%  through  the 
wall  thickness  in  the  100%  overstrain  specimens. 

The  effect  of  large  plastic  zones  on  crack  growth  rate 
in  aluminum  has  been  studied  by  Elber^.  He  has  found 
that  crack  closure  can  occur  in  a zero  to  tensile  fatigue 
cycle  before  all  the  tensile  load  is  removed , due  to  the 
compressive  residual  stresses  in  the  plastic  zone.  This 
means  that  the  effective  crack  opening  fatigue  cycle  is  de- 
creased; hence,  a slower  crack  growth  rate.  Based  on  this 
closure  model,  the  conclusion  can  be  drawn  that  large  plastic 
zones  retard  crack  growth,  yet  in  the  above  experiment,  it 
seems  that  large  plastic  zones  accelerate  crack  growth. 

A similar  model  can  be  used  to  explain  the 
accelerated  crack  growth  that  occurred  in  the  disk  experi- 
ment. The  tensile  residual  stress  causes  a crack  under  no 
applied  load  to  be  open,  with  stress  intensity  factor  of: 

^^min  °R  (18) 

After  application  of  the  load,  the  stress  intensity 
factor  is: 

(Kl)max  a°R  + °APP  H9) 
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If  the  sum  of  the  residual  and  applied  stresses  exceeds 
the  yield  strength,  the  stress  intensity  factor  at  maximum 
load  is: 

(Kl)maxa°y  (20) 

When  the  load  is  removed  after  yielding  has  occurred, 
some  of  the  tensile  residual  stress  is  locally  relaxed  by 
an  amount  0*.  which  causes  the  crack  to  close  more  under 
minimum  load.  The  stress  intensity  factor  at  minimum  load 
i s : 

(KI>min  “ °R  - o*  (21) 

When  the  load  is  again  applied,  the  stress  intensity 
factor  at  maximum  load  is: 

(Kl)max  a°R  - ^ ♦ °APp  (22) 

If  the  sum  of  applied  and  residual  stresses  still  ex- 
ceed the  yield  strength,  the  stress  intensity  factor  at 
maximum  load  again  is: 

(Kt) 

max  “ °y  (23) 

Since  the  stress  intensity  factor  at  maximum  load  is  con- 
stant (equations  20,23),  but  at  minimum  load  is  reduced  pro- 
portional to  o*  (equation  22),  the  range  of  stress  intensity 
factor  (A  K i ) is  increased,  causing  faster  crack  growth  rate. 
If  a*  is  large  enough,  a large  increase  in  crack  growth 


rate  is  expected. 


Since  the  plastic  zone  appeared  at  shallower  crack 
depths,  and  was  larger  in  the  100%  overstrain  specimens 
than  the  50%  overstrain  specimens,  o*  associated  with  100% 
overstrain  is  expected  to  be  larger  than  a*  associated 
with  50%  overstrain;  therefore,  faster  crack  growth  rate 
is  expected  with  more  autofrettage. 

VI . CONCLUSION 

Autofrettage  has  a deleterious  effect  on  fatigue  crack 
propagation  in  externally  flawed  thick-walled  cylinders. 

The  effect  is  believed  to  be  the  result  of  a crack  closure 
phenomenon  caused  by  large  plastic  zones  accompanying  crack 
growth  which  effectively  increases  the  stress  intensity 
factor  range.  If  the  plastic  zones  could  be  reduced,  by 
experimenting  with  thicker  specimens,  the  large  increase  in 
crack  growth  rate  would  probably  also  be  reduced. 

VII.  ACKNOWLEDGEMENT 

The  author  would  like  to  acknowledge  the  help  of  Mr.  W. 
Mortimer  for  his  assistance  in  conducting  the  above  experiments, 
Mr.  J.  Bak  for  obtaining  the  specimens.  Dr.  M.  Hussain,  Mr. 

W.  Lorenson  and  Mr.  G.  Pflegl,  for  their  suggestions  and 
assistance  in  the  use  of  the  program  NASTRAN,  and  Mr.  J.  Under- 
wood for  his  many  suggestions,  guidance  and  assistance  in 
analyzing  the  data,  and  Mrs.  M.  Dietz  for  typing  the  manuscript. 


38 


REFERENCES  : 


1.  Timoshenko,  S.,  "Strength  of  Materials,  Part  II",  Third 
Edition,  D.  Van  Nostrand  Co.,  Princeton,  N.J.,  1956. 

2.  Davidson,  T.E.  and  Kendall,  D.P.,  "The  Design  of  Pressure 
Vessels  for  Very  High  Pressure  Operation",  Watervliet 
Arsenal  Technical  Report,  WVT-6917,  1969. 

3.  Paris,  P.C.,  "The  Fracture  Mechanics  Approach  to  Fatigue" 
in  "Fatigue,  An  Interdisciplinary  Approach",  Burke,  J.J., 

Reed,  N.L.,  Weiss,  V.,  Editors,  Syracuse  University  Press, 

(1964),  pp  107-132. 

4.  Kobayashi,  A.S.,  "A  Simple  Procedure  for  Estimating  Stress 
Intensity  Factor  in  Region  of  High  Stress  Gradient", 

Interim  Technical  Report  No.  1,  U.S.  Army  Research  Grant 
No.  DA- ARO- D- 3 1 - 1 24 - 73- 638  , 1973. 

5.  Chopra,  P.S.,  "Finite  Element  Fracture  Mechanics  Analysis 
of  Creep  Rupture  of  Fuel  Element  Cladding",  Nuclear 
Engineering  and  Design  29,  1974. 

6.  Kobayashi,  A.S.,  Maiden,  D.E.,  and  Simon,  B.J.,  "Application 
of  Finite  Element  Analysis  Method  to  Two-Dimensional 
Problems  in  Fracture  Mechanics",  ASME  Paper  No.  69-WA/PVP- 1 2 , 

1969,  pp  1-7. 

7.  Emery,  A.F.  and  Segedin,  "The  Evaluation  of  the  Stress 
Intensity  Factors  for  Cracks  Subjected  to  Tension,  Torsion, 
and  Flexure  by  an  Efficient  Numerical  Technique",  Journal 

of  Basic  Engineering,  pp  387-393,  June  1972.  I 

8.  Wigglesworth , L.A.,  "Stress  Distribution  in  a Notched  Plate", 
Mathematica,  4 (157),  76-79. 

9.  Shannon,  R.W.E.,  "Crack  Growth  Monitoring  by  Strain  Sensing", 

Pres.  Ves.  and  Piping  (1),  pp  61-73,  1973.  j 

10.  Hussain,  M.A.,  Lorensen,  W.E.,  and  Pflegl,  G. , "The  Quarter- 
Point  Quadratic  Isoparametric  Element  as  a Singular  Element 
for  Crack  Problems",  to  be  published  in  NASA  Technical 
Memorandum,  NASTRAN : Users  Experiences,  October  1976. 


39 


• ail* 


» .»• 


■1  k . Jl 


11.  Wes teagaard , H.M.,  "Bearing  Pressures  and  Cracks", 
Transactions  of  the  ASME,  Journal  of  Applied  Mechanics, 

1939  . 

12.  Underwood,  J.H.,  private  communication,  July  1976. 

13.  Hussain,  M.A.,  Pu,  S.L.,  and  Underwood,  J.,  "Strain  Energy 
Release  Rate  for  a Crack  Under  Combined  Mode  I and  Mode  II", 
Fracture  Analysis,  A STM  STP560,  American  Society  for  Testing 
and  Materials,  pp  2-28,  1974. 

14.  Clark,  W.G.,  "How  Fatigue  Crack  Initiation  and  Growth 
Properties  Affect  Material  Selection  and  Design  Criteria", 
Metals  Engineering  Quarterly,  pp  16-22,  August  1974. 

15.  Elber,  Wolf,  "The  Significance  of  Fatigue  Crack  Closure", 
Damage  Tolerance  in  Aircraft  Structures,  ASTM  STP486  , 
American  Society  for  Testing  and  Materials,  pp  230-242, 

1971  . 


40 


WATERVLIET  ARSENAL  INTERNAL  DISTRIBUTION  LIST 


May  1976 
No.  of  Copies 


COMMANDER  1 

DIRECTOR,  BENET  WEAPONS  LABORATORY  1 

DIRECTOR,  DEVELOPMENT  ENGINEERING  DIRECTORATE  1 

ATTN:  RD-AT  1 

RD-MR  1 

RD-PE  1 

RD-RM  1 

RD-SE  1 

RD-SP  1 

DIRECTOR,  ENGINEERING  SUPPORT  DIRECTORATE  1 


2 

1 
1 
1 
1 

TECHNICAL  LIBRARY  5 

TECHNICAL  PUBLICATIONS  6 EDITING  BRANCH  2 


DIRECTOR,  RESEARCH  DIRECTORATE 
ATTN:  RR-AM 

RR-C 
RR-ME 
RR-PS 


DIRECTOR,  OPERATIONS  DIRECTORATE  1 
DIRECTOR,  PROCUREMENT  DIRECTORATE  1 
DIRECTOR,  PRODUCT  ASSURANCE  DIRECTORATE  1 
PATENT  ADVISORS  1 


EXTERNAL  DISTRIBUTION  LIST 


OFC  OF  THE  DIR.  OF  DEFENSE  RSE 
ATTN:  ASST  DIRECTOR  MATERIALS 

THE  PENTAGON 
WASHINGTON,  D.C.  20 3 IS 

CDR 

IJS  ARMY  TANK-AUTMV  COMD 
ATTN:  AMDTA-UL 

AMSTA-RKM  MAT  LAB 
WARREN,  MICHIGAN  48090 

CDR 

P1CATINNY  ARSENAL 
ATTN:  SARPA-TS-S 

SARPA-VP3  (PLASTICS 
TECH  EVAL  CEN) 

DOVER,  N.J  07801 


CDR 

US  ARMY  ARM!  COMD 
ATTN:  AMSAR-PPW-IR 

AMSAR-RD 
AMSAR-RDG 

ROCK  ISLAND,  IL  61201 
CDR 

US  ARMY  ARMT  COMD 
FLD  SVC  DIV 
ARMCOM  ARMT  SYS  OFC 
ATTN:  AMSAR-ASF 
ROCK  ISLAND.  IL  61201 

CDR 

US  ARMY  ELCT  COMD 
FT  MONMOUTH,  NJ  07703 


CDR 

FRANKFORD  ARSENAL 
ATTN : SARFA 

PHILADELPHIA,  PA  19137 

DIRECTOR 

US  ARMY  BALLISTIC  RSCH  LABS 
ATTN:  AMXBR-LB 

ABERDEEN  PROVING  GROUND 
MARYLAND  21005 

CDR 

US  ARMY  RSCH  OFC  (DURHAM) 
BOX  CM,  DUKE  STATION 
ATTN:  RDRD-IPL 

DURHAM,  NC  27706 


CDR 

REDSTONE  ARSENAL 
ATTN:  AMSMI-RRS 

AMSMI-RSM 
ALABAMA  35809 

CDR 

ROCK  ISLAND  ARSENAL 
ATTN:  SARRI-RDD 

ROCK  ISLAND,  IL  61202 

CDR 

US  ARMY  FGN  SCIENCE  5 TECH  CEN 
ATTN : AMXST-SD 

220  7TH  STREET  N.E. 
CHARLOTTESVILLE,  VA  22901 


CDR  DIRECTOR 

WEST  POINT  MIL  ACADEMY  US  ARMY  PDN  EQ.  AGENCY 

ATTN:  CHKN,  MECH  ENGR  DEPT  ATTN:  AMXPE-MT 

WEST  POINT,  NY  10996  ROCK  ISLAND,  IL  61201 


CDR 

HQ,  US  ARMY  AVN  SCH 
ATTN:  OFC  OF  THE  LIBRARIAN 

FT  RUCKER,  ALABAMA  36362 


1 copy  to  each 


EXTERNAL  DISTRIBUTION  LIST  (Cont) 
2 copies  to  each 


OR 

IIS  NAVAL  KPNS  LAB 
CHIEF,  MAT  SCIENCE  DIV 
ATTN:  MR.  I).  MALYEVAC 

hahlgren.  va  22448 
mui  ci  or 

NAVAL  RSCH  LAB 
ATTN:  (HR.  MECH  DIV 
W ASH  1 NCTON , l).C.  2037S 

1>I  RECTOR 
NAVAL  RSCH  LAB 

coni:  2d-27  (i)ocu  lib.) 

WASHINGTON,  D.C.  20375 

NASA  SCIENTIFIC  5 TECH  INFO  FAC 
PO  BOX  8757,  ATTN:  ACQ  BR 

BALTIMORE/WASHINGTON  INTL  AIRPORT 
MARYLAND  21240 


CDR 

US  ARMY  MOB  EQUIP  RSCH  6 DEV  COMD 
ATTN:  TECH  DOCU  CEN 
FT  BELVOIR,  VA  22060 

CDR 

US  ARMY  MAT  RSCH  AGCY 
ATTN:  AMXMR  - TECH  INFO  CEN 

WATERTOWN,  MASS  02172 

CDR 

WRIGHT- PATTERSON  AFB 
ATTN : AFML/MXA 

OHIO  45433 

CDR 

REDSTONE  ARSENAL 

ATTN:  DOCU  6 TECH  INFO  BR 

ALABAMA  35809 


DEFENSE  METALS  INFO  CEN 
BATTELLE  INSTITUTE 
505  KING  AVE 
COLUMBUS,  OHIO  43201 

MANUEL  E.  PRADO  / G.  ST1SSER 
LAWRENCE  LIVERMORE  LAB 
PO  BOX  808 
LIVERMORE,  CA  94550 

l)R.  ROBERT  QUATTRONE 
CHIEF,  MAT  BR 
US  ARMY  R&S  GROUP,  EUR 
BOX  65,  FPO  N.Y.  09510 


12  copies 
CDR 

DEFENSE  DOCU  CEN 
ATTN:  DDC-TCA 

CAMERON  STATION 
ALEXANDRIA,  VA  22314 


NOTE:  PLEASE  NOTIFY  CDR,  WATERVLIET  ARSENAL,  ATTN:  SARWV-RT-TP, 

WATERVLIET,  N.Y.  12189,  IF  ANY  CHANGE  IS  REQUIRED  TO  THE  ABOVE. 


. /■'  ‘ ■rV  ? 


■J  AaVil.-  -v- 


